Luttinger Liquid Physics in the Superconductor Vortex Core 
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We study several aspects of the structure of vortices in conventional s-wave Type II supercon- 
ductors. It is well-known that there are low energy quasiparticles bound to the core of a vortex. 
We show that under certain conditions, these quasiparticles form a degenerate Fermi gas with a 
finite density of states at the Fermi energy. In three dimensional superconductors, the result is 
a one dimensional Fermi gas of quasiparticles bound to the core of a vortex line. As is usual in 
one dimensional systems, interactions between the quasiparticles lead to Luttinger liquid behaviour. 
This may be probed through STM tunneling into the vortex core. We further suggest that a novel 
Peierls-type instability in the shape of the vortex line may develop due to interactions between the 
quasiparticle gas and fluctuations in the vortex line shape. 



I. INTRODUCTION 



The electronic properties of Type // superconductors 
in the mixed phase are known to have several interesting 
features. In some pioneering work, Caroli et. al. pre- 
dicted theoretically the existence of low energy quasipar- 
ticle states bound to the core of a single isolated vortex in 
an s-wave superconductor. These low lying states are free 
to move parallel to the vortex, and so give rise to a series 
of one dimensional bands (minibands) as shown in Fig 
1(a). The system is still gapped, though the energy gap 
to the lowest excitation (the minigap) is typically much 
smaller than the bulk gap. Subsequent theoretical work 
[^-Q helped to put this prediction of low energy quasi- 
particle states bound to the vortex core on firm ground. 
Striking experimental evidence for the existence of these 
bound states of the vortex was provided by the scanning 
tunneling microscopy (STM) studies of Hess et al [|. 
These experiments were able to image the vortex quasi- 
particle states as a function of their energy, and agreed 
well with theory. 

Most theoretical work on these low energy quasiparti- 
cle states in the vortex core has neglected the effect of 
the Zeeman coupling between the magnetic field and the 
spin of the quasiparticles. However as discussed in Ref. 
Q , and as we argue below, the Zeeman coupling plays an 
important role in determining the low temperature prop- 
erties of the system. The Zeeman coupling splits the 
quasiparticle energy levels, and the minigap is decreased 
since the energy of one spin species is lowered towards 
the Fermi energy. For sufficiently strong Zeeman split- 
ting, the minibands of one spin species will be brought 
down below the Fermi energy, and a filled Fermi sea of 
spin polarised quasiparticles is formed. We show that the 
magnetic fields needed in typical materials to form this 
degenerate quasiparticle system are not large, and could 
be much smaller than Hc2- 

We consider the effect of quasiparticle-quasiparticle in- 
teractions (ignored in the BCS mean field theory) on the 



low energy properties of the quasiparticles in the core 
of the vortex. The quasiparticles are bound in the di- 
rection perpendicular to the vortex line but are free to 
move along it, thus providing an interesting realization of 
a one dimensional system inside the superconductor. It 
is well-known in the theory of normal metals that inter- 
action effects are dramatic in one dimension: the generic 
ground state of the interacting electron system is not a 
Fermi liquid but a different beast, the Luttinger liquid. 
We therefore focus attention primarily in the regime of 
magnetic fields well below Hc2 where the vortices may be 
treated in isolation. Are some of the striking properties 
of interacting ID Fermi systems, such as Luttinger liquid 
physics, also present in the vortex quasiparticle system? 
The answer is no, if the interactions and the Zeeman 
coupling are weak. The presence of the (mini)gap im- 
plies that the T=0 state of this system is quite insensi- 
tive to weak interactions. However if the Zeeman energy 
is large enough to start filling the miniband, the system 
is gapless - interaction effects are then crucial, and we 
argue that the system is correctly described as a spin 
polarised Luttinger liquid. We then find that the in- 
teraction strength g which controls the Luttinger liquid 
exponents is a function of the miniband filling, and in 
principle could even be negative. Luttinger liquids have 
several interesting properties - the fermion correlation 
functions are power law with anomalous exponents and 
the single particle density of states at the Fermi points 
vanishes as a power law. There have been a small num- 
ber of experimental realizations of Luttinger liquid be- 
haviour in systems such as ID semiconductor wires 
fractional quantum hall effect edge states |^ and carbon 
nanotubes ||l^. We show that a Luttinger liquid can be 
realised under appropriate conditions in the vortex core. 
Tunneling into a Luttinger liquid leads to a character- 
istic power law tunneling conductance, and we propose 
an experiment involving Scanning Tunneling Microscopy 
(STM) measurements on the vortex core as a probe of 
Luttinger liquid behaviour of the vortex quasiparticles. 
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The quasiparticles also interact with the collective 
modes of the vortex, which can have interesting conse- 
quences. In particular, if inter-vortex interactions (which 
become increasingly important as we increase the field 
towards i?c2) are sufficiently strong, we find that a vor- 
tex analog of the Peierls effect could occur. The vortex 
line spontaneously undergoes a periodic modulation of 
its profile to lower the energy by opening up a gap in the 
quasiparticle spectrum. 

The rest of this paper is organised as follows. In Sec- 
tion ||, we discuss in a little more detail the core states 
of an isolated vortex, and the formation of the degen- 
erate quasiparticle g as i n the presence of the Zeeman 
splitting. In Section III we show that interactions be- 



tween the quasiparticles can lead to them forming a Lut- 
tinger liquid, with a varying Luttinger liquid exponent 
and consider how this state may be observed in STM ex- 
periments. In Section IV, we couple the quasiparticles to 
the vortex collective modes and discuss in particular the 
possibility of a vortex Peierls transition when inter- vortex 
interactions are present. Section discusses the validity 
of some of our approximations and highlights directions 
where more work needs to be done. 



II. VORTEX QUASIPARTICLE STATES IN THE 
NON-INTERACTING LIMIT 

In this section we review the spectrum of the vortex 
core states and discuss the effect of the Zeeman coupling 
on it. We begin with the BCS model Hamiltonian in the 
presence of an arbitrary magnetic field B — V x A and 
pair potential A(r) (wc will include the Zeeman term 
later): 

Hbcs = J d'^x 4(a;)e[-iV - eA{x)]ca{x) + 

4 {x)^{x)c\ {x) + q (a;) A* {x)c^ {x) (1) 

where e[p\ is the kinetic energy measured from the chem- 
ical potential. For simplicity, we will assume a quadratic 
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energy dispersion e{p) = ^ — Ep- It is convenient to 
rewrite the Hamiltonian using the change of variables 

di{x) = c^{x) 

dl{x) = c:^{x) 

When written in these variables, there are no anoma- 
lous terms in the Hamiltonian: 



H 



BCS 



d\{S,{x) 4(^))^( dilx) 



H = 



(2) 
(3) 



e[-iV-eA{x)] A*(a;) 

A(a;) -e[iV - eA{x) 

Physically, the number of d-particles corresponds to 
the z-component of the electron spin. The absence of 
anomalous terms in the d representation thus reflects con- 
servation of the z-component of the electron spin in the 
BCS Hamiltonian of Eqn. 0. 



The Hamiltonian is diagonalised by going over to new 
quasiparticle operators that are defined as: 



7L = j {ua{x)d\{x) + Va{x)dl{x))d'^x 
7L = / {v:ix)d\ix)-ulix)dUx))d''x 



The functions Ua{x) and Va(x) are found by solving the 
following eigenvalue equation (with Ea > 0): 



Vaix) J "I Vaix) 



(4) 



This is just the Bogoliubov dcGcnncs (BdG) equation 
for the quasiparticle states. The BdG Hamiltonian can 
be recast in a compact matrix notation as: 

H ^ e[-iVT^ - eA{x)]T^ + A*{x)t+ + A{x)t- 

where r are the 2x2 Pauli matrices. In terms of the 
quasiparticle operators 7, the Hamiltonian is simply: 



Hbcs = ( "fLi^) iLix) 



Ea \ f lla{x) 



~Ea 



l2a{x) 



(5) 



The ground state has all of the negative energy states 
filled so 72q|0) = 0; all positive energy states are unoccu- 
pied 7iq|0) — 0. Note that a spin t excitation is created 
by (7i ) while a spin | excitation is created by (72 ) acting 
on the ground state. 

After these generalities, let us specialise to the quasi- 
particle states bound to the vortex line in a superconduc- 
tor. Consider an isolated, straight vortex oriented along 
the z-axis, sitting in a clean, conventional Type II super- 
conductor. This is realised if the applied magnetic field 
is just above Hd- The vortex bound state energy levels 
are obtained by solving Eqn. ^ with the appropriate gap 
profile (A(r) = |A(r)|e*'^, where (p is the azimuthal angle 
about the vortex). It is useful to perform the singular 
gauge transformation that makes A real everywhere (the 
London gauge): 



(6) 
(7) 



H' = [^(-*Vr^ -eA- ^ f - E^jr^ + |A(r)|r- (8) 



The transformed quasiparticle wavefunctions 



obey antiperiodic boundary conditions to keep ( ) sin- 



gle valued. 



E 



Fermi sea of quasiparticles 



3/2 5/2 




FIG. 1. Energy level structure of the states bound to the core of a vortex (a) Ignoring the effect of the magnetic field. The 
negative energy states are occupied and shown with the thick line. Energy is measured in units of the minigap energy (« ■^) 
and (b) In the presence of a sufficiently strong magnetic field. The Zeeman splitting causes the = \ miniband to start being 
occupied and a degenerate Fermi sea is formed. 



We first review the results of Rcf. on the struc- 
ture of an isolated vortex. In Ref. [|l| the effect of the 
magnetic field on the low energy quasiparticles is not in- 
cluded. This is only a good approximation for an isolated 
vortex in an extreme Type II superconductor; then the 
flux enclosed in the spatial region where the low energy 
states are bound is a small fraction of the flux quantum, 
and hence the magnetic field has little effect on these 
states. Subsequently we will see how these results need 
to be modified on including the effect of the magnetic 
field. The low energy spectrum of the vortex core states 
then is : 



k 771 



(9) 



and is shown in Fig 1(a). The corresponding quasiparti- 
cle wavefunctions in cylindrical coordinates take the form 



9,1k At) 



(10) 



where ji G |, §...} appears as an angular momentum 
in the wavefunctions and is a half integer due to the anti- 
periodic boundary conditions, is the momentum along 
the vortex line, c is a constant of order one and ^ is the 
zero temperature coherence length. We have denoted by 
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Aoo the value of the gap very far away from the vortex. 
The kz dependence of the energy comes from the motion 
along the vortex line. The energy scale of the minigap is 

-g^ and is much smaller than the bulk gap Aqo ■ The ra- 
dial wave funcions / and g are mainly confined to within 
a radius ^ (the coherence length) about the vortex for 
the low energy states fc^) ^ Aqo)- 

Now consider the effect of the magnetic field B — Bz. 
We assume though that the field is much less than Hc2 
so that the vortices are well separated. Two effect oc- 
cur Q : (a) the Zeeman term splits the energy levels and 
(b) there is an orbital effect. We consider these in turn, 
starting with the orbital effect. 

For the magnetic fields of interest, the inter- vortex sep- 
aration is typically smaller than the penetration depth 
so that the field in the superconductor is fairly uniform 
and equal to the external field {— Bz, and we take 
A{x) = y, x, 0)). In the presence of the vector po- 
tential additional terms are generated in the BdG Hamil- 
tonian which in the London gauge take the form: 



5H = 5Hi + 5H2 
Zm 



jeAjx))' , 



\e\B^ 
2m ' 



(11) 
(12) 

(13) 



5H2 may be neglected since the first term is small com- 
pared to the minigap for fields well below Hc2 and the 
second term only makes a small shift of the Fermi energy. 
So we are left with 6H1 which can be written as: 



e e 
Zm Zm 



(14) 



This does not affect the eigenstates (|T^) but shifts their 
energy by 



E{fi,kz,B)^E°{fi,K) + {^- 



e\B 



2m 



(15) 



The orbital effect therefore increases the magnitude of 
energy of the quasiparticle states. 
The Zeeman term takes the form 



Hz = 



2 



ic\ix)c^{x)^c\{x)c^ix))d''x (16) 
dl{x)di{x) + dl{x)d2{x)) d'^x (17) 



(18) 



(In the second line, we have dropped an irrelevant con- 
stant and fJ,B = is the Bohr magneton). The Zee- 
man term thus behaves as a 'chemical potential' for the 
d-particles Increasing the magnetic field raises this 
chemical potential. Beyond a certain field it enters the 



first miniband and a degenerate Fermi sea of quasiparti- 
cles forms as shown in Fig 1(b). The condition for this 
is simply: 



Ei^=l,kz=0,B)<iI^ 



or equivalently 



E^i^,^^,kz^o)<^-^^^,BB 



(19) 



(20) 



Thus a degenerate gas of spin polarised fermions (in 
this case all spin t) can form due to the Zeeman coupling 
if the magnetic field exceeds Hcz = ^-E'°(m = = 0) 
(assuming g = 2). Note that it is only the first miniband 
that can be brought below the chemical potential. For 
/i > i, the orbital effect wins over the Zeeman splitting, 
and the levels do not cross the chemical potential. 

The magnetic fields needed to begin filling the mini- 
band states are very reasonable, since the minigap energy 
is typically small. To illustrate this point we consider the 
case of NbSe2 for which numerical calulations of the vor- 
tex quasiparticle spectrum are available Q. NbSe2 has 
a superconducting transition temperature of Tc — 7.2K, 
and numerical calculations find the minigap to be 0.4 
Kelvin. The magnetic field required to close the mini- 
gap is [Hcz] is 1-3 Tesla (we have assumed g = 2) while 
Hc2 is larger at 3.2T for this material. Better candidate 
materials will have smaller ratios of Uf^ than NbSco. 

In passing we note that the Zeeman splitting gives rise 
to quasiparticles even at zero temperature, and in doing 
so 'melts' some of the condensate at the centre of the 
vortex. The new profile A(r), can be calculated from 
the self consistency equation. This is analogous to the 
Kramer-Pesch effect ||^, where temperature does the job 
of melting the condensate and expanding the vortex. We 
shall not take into account the effects of the changing gap 
profile on the quasiparticle states, as they are expected 
to be small. 



III. INTERACTION EFFECTS AND LUTTINGER 
LIQUID FORMATION 

Interactions between fermions in clean ID systems 
have a dramatic effect resulting in Luttinger liquid be- 
haviour. Here we shall consider the effect of interactions 
between the quasiparticles in the vortex core that form 
the ID Fermi gas. A few comments on the nature of 
these interactions [ p^ is in order. The underlying mi- 
croscopic electronic system is characterized, at low en- 
ergies, by three qualitatively different kinds of interac- 
tions: (a) the spin singlet density-density interaction in 
the particle-hole channel (b) the triplet spin density-spin 
density interaction in the particle-hole channel and (c) 
the interaction in the particle-particle channel. Conven- 
tional BCS superconductors arise when the interactions 
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in the particle-particle channel are attractive. Indeed the 
electron gas is unstable toward the formation of Cooper 
pairs in the presence of such an attractive interaction. 
The BCS theory simply treats this attractive interaction 
in the particle-particle channel in a mean field approx- 
imation. One still however has to deal with the inter- 
actions in the particle-hole channel. (These are simply 
ignored in the so-called reduced BCS Hamiltonian). In- 
clusion of these leads to interactions between the quasi- 
particles in the superconductor. Note that the long- 
ranged Coulomb repulsion between the electrons in the 
singlet channel is screened out by the condensate in the 
superconducting phase. Thus the residual interactions 
between the quasiparticles are short-ranged. A further 
source of interactions between the quasiparticles comes 
from the inclusion of fluctuations of the mean field order 
parameter (again usually ignored in BCS theory). 

Our general conclusions are insensitive to the precise 
form of this short-ranged interaction between the quasi- 
particles. We will therefore, for illustrative purposes, 
consider a particular model interaction. We add to the 
BCS Hamiltonian Eqn.|l] the "interaction" term 

H^nt^ J ^F(|x-a;'l)ct(x)ct,(x')c.'(a^')c.(a^)rf'^c?^^' 

era' 

(21) 

where V{\x—x']) is assumed to be short-ranged. To study 
the effect on the vortex core states, it is convenient to re- 
express this in terms of the 7 operators introduced in 
Section Q 




x{vp'{x')jip' - Up,{x')-f2p') ■■ d'^xd'^x' 
+ ... + ... + ... 

where the normal ordering is with respect to the super- 
conductor vacuum, and we have written out only the 
first of the four terms present from different spin combi- 
nations. 

Let us study the effect of interactions when we may 
treat the vortices as isolated and the lowest miniband is 
occupied with the two Fermi points at ±fcy (so q > 
Ezeeman > El q). Noticc that the Fermi Sea consists 
of spin up quasiparticles only. The excitation spectrum 
for this case contains a low energy part of 'particle-hole' 
excitations in the vicinity of the Fermi points, and a high 
energy sector that involves either 71 quasiparticles hop- 
ping into higher minigap states, or the destruction of 72 
quasiparticles from deep below the chemical potential. 
The latter two are gapped, with the gap being of order 
the minigap energy (« jpr). Since we are interested in 



the effects of weak interactions on the T = state of 
the system, we confine ourselves to the low energy sec- 
tor of the problem. Formally, we can achieve this by 
defining a projection operator that only retains states 
close to the Fermi points, i.e. {71 quasiparticles, /i = ^, 
kz e (fc/ - e,kf + e) or £ (— fc/ - e, -fc/ + e); e < fc/}. 
Now, we can project the interaction term Eqn(p^) down 
into this subspace, and it is easy to see that the only 
non-trivial term that remains is: 

Hf„,^uJ2pLiq)PR{~q) (23) 

where the left (and similarly the right) density opera- 
tors PL and PR are constructed from the fermion opera- 
tors in the usual way, pL for example is given by: 

\k+kf\<e 

and Ffc = Ji 1 i^. In general, doing a quantitative cal- 
culation of U for a real material is not an easy task be- 
cause a complete knowledge of the interaction potential 
and quasiparticle wavefunctions is needed. We simply 
note that U will be a function of the filling of the mini- 
band (U = U{kf )), and in principle can even be negative 
(attractive interaction between the quasiparticles) p3[ |. 

For the particular interaction of Eqn. ^ we can ex- 
press U in terms of the interaction potential and the 
wavefunction of the states near the Fermi points: 

f V{x-x\y-y',z)Qk,{x,v)Qk,{x',v')[l-e^'''' 

(24) 

Here the quasiparticle 'charge' Qkf is defined by 

Qkfix,y) = \ui j^^{x,y)\'^ - |f i,fc^(x,y)p (25) 

We thus have a system of one-dimensional fermions 
with two Fermi points interacting through H^^. The 
fermions are spin-polarized (and hence, effectively, spin- 
less). It is well-known that such a system is correctly de- 
scribed not as a Fermi liquid, but rather as a Luttinger 
liquid. 

Luttinger liquids are characterised by power law cor- 
relations. The exponents are controlled by the dimen- 
sionless quantity g — where f/ is the Fermi ve- 

locity of the one dimensional system. For example, 
(r'^(z)r(z')) « {l/\z — z'Y) where z is the position along 
the vortex line and the exponent v — (Xj \J\ — g^). Re- 
latedly, tunneling into the bulk of a Luttinger liquid is 
characterised by nonlinear I-V characteristics, in fact at 
low voltages the tunneling current-voltage relation is of 
the form / oc . This suggests that STM could be used 
to experimentally verify Luttinger liquid formation in the 
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vortex cores under the conditions described above. At 
sufficiently low temperatures when the Luttinger liquid 
behaviour is dominant, we expect the tunneling conduc- 
tance close to the centre of the vortex to be given by 

(26) 

and Vf, is the exponent for tunneling into the edge of a 
Luttinger liquid i^e — (1+3)^(1—3)"^ [0- The temper- 
atures at which such Luttinger liquid behaviour will be 
observed is necessarily small compared to the degeneracy 
temperature, which for the half filled miniband of NhSei 
is ~ 0.4Kelvin These effects may be observed at 

higher temperatures by using a superconductor with a 
larger minigap, but this would also require larger mag- 
netic fields to close the minigap via the Zeeman splitting. 

IV. COUPLING TO VORTEX COLLECTIVE 
MODES AND VORTEX PEIERLS EFFECT 

So far we have only considered straight and rigid vor- 
tex lines, however in reality the vortex is a soft object and 
can undergo shape fluctuations. These collective modes 
of the vortex interact with the quasiparticles bound to 
the vortex line. As we show below the interacting quasi- 
particle - collective mode system is similar to the prob- 
lem of interacting electrons and phonons in one dimen- 
sion, which is known to undergo a Peierls transition (for 
commensurate filling) in which the lattice spontaneously 
distorts and a gap opens at the Fermi points. The gain in 
electronic energy by opening of the gap offsets the elastic 
energy cost of the lattice distortion. The electron scat- 
tering from the vicinity of one Fermi point to the other 
(2fci? scattering) involves vanishing energy denominators 
and is responsible for this instability. 

We can ask if an analogous effect occurs in the degener- 
ate ID quasiparticle system in the vortex core interacting 
with the shape fluctuations of the vortex - which can give 
rise to 2kp scatterings [|7j. Note that for an isolated 
vortex line, which is a one dimensional object it is not 
possible to spontaneously break the continuous symme- 
try of translations along the vortex line. However in the 
vortex lattice if inter vortex interactions are sufficiently 
strong, a vortex analog of the Peierls effect could occur 
at finite temperature. This issue was raised earlier by 
Bouchaud who modelled the vortex core simply as a 
wire of normal electrons. However we believe this is not 
an adequate model to discuss this phenomena. A degen- 
erate Fermi system is formed in the vortex only in the 
presence of the Zeeman splitting. 

The collective modes of the vortex may be classified ac- 
cording to their angular momentum. Scalar modes scat- 
ter quasiparticles within the same [i miniband, while the 
lowest order coupling of quasiparticles to collective modes 



of higher angular momenta (say X) will scatter quasipar- 
ticles from one miniband {\x) to a different one (yit ± /) 
by conservation of angular momentum. Since the low 
energy excitations all lie in a single [i = \ miniband we 
only consider coupling of quasiparticles to the scalar col- 
lective modes of the vortex. 

Scalar vortex modes modulate the radial profile of the 
vortex as one moves along the vortex line. For instance, 
a periodic modulation of the vortex radius (defined say 
as when the order parameter reaches half its asymptotic 
value) is an example of a scalar deformation of the vor- 
tex. If the scalar normal modes of the vortex are labelled 
by n, let <Jn{qz) create the n*'' radial mode with momen- 
tum along the vortex line. The coupling of the vortex 
modes to the quasiparticles (of the partially filled mini- 
band) can be derived from the BdG equation, or simply 
from symmetry is found to be: 

H!:?'""- = E 9^ntik, q, n)a„(g)r^^r, + h.c. (27) 

q,k,n 

The Hamiltonian for the scalar modes is: 

Hmodes = ^ E{n, q)(j\{q)an{q) (28) 

Thus, scalar modes of the vortex of momentum 2k f will 
scatter quasiparticles from one Fermi point to the other. 
As we have noted previously, for an isolated vortex this 
does not lead to a gapped phase and the Luttinger liquid 
behaviour will persist. However if inter vortex interac- 
tions are sufficiently strong, then a finite temperature 
transition to a vortex Peierls phase can occur. Which of 
the modes has the highest transition temperature is 
a function of gint^ energy and inter vortex coupling for 
that mode. In the vortex Peierls phase, the quasiparti- 
cle spectrum is gapped, and we have a collective mode 
condensate i.e. (cr„o(2fc/)) ^ which implies a periodic 
modulation of the vortex profile. 

We do not attempt to predict the precise conditions for 
the experimental realization of this phase. Increasing the 
vortex density, and hence increasing intervortex interac- 
tions can raise the transition temperature, as long as the 
core quasiparticles do not hop between vortices and de- 
stroy the one-dimensional nature of the system. Within 
a mean field approximation, the Fermi gas always under- 
goes the Peierls transition. The instability is enhanced 
for repulsive Luttinger liquids ]l9t . Thermal transport 
measurements along the vortex line are expected to be 
sensitive to the vortex Peierls transition. 

The Luttinger liquid behaviour of the vortex quasipar- 
ticles would also be destroyed if pairing of quasiparticles 
were to occur. This would induce an additional (eg. p- 
wave) pairing amplitude in the vortex core, and could 
also lead to breaking of the rotational symmetry of the 
vortex line. This interesting situation was considered by 
Makhlin and Volovik in [pO|. Once again, the transition 
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temperature for this instability is nonzero only in the 
presence of suitable inter-vortex interactions. 



V. DISCUSSION 

We have shown that a one-dimensional degenerate gas 
of quasiparticles can form in the vortex core, as a result 
of the Zeeman coupling and can serve as a laboratory for 
one dimensional physics. We showed that quasiparticle 
interactions drive the formation of a Luttinger liquid in- 
side the vortex core. The interaction strength U of the 
Luttinger liquid was shown to be a function of the filling, 
and could even be negative. Thus we have a ID inter- 
acting Fermi system, where novel features arise due to 
the fact that the fermions involved are superconducting 
quasiparticles and the 'wire' confining them is a vortex. 

We now examine in more detail some of the approxi- 
mations that have been made, and the prospects for ver- 
ifying the results we predict in experiments. 

One of the main approximations that we have made 
is to treat each vortex as effectively isolated, which gives 
rise to the one dimensional nature of the vortex quasipar- 
ticle system. As the magnetic field is increased towards 
Hc2, the vortices get closer to each other and the wave- 
fuctions of the vortex core states start to overlap. This 
leads to quasiparticle hopping between vortices that will 
eventually destroy the one dimensional nature of the sys- 
tem. To keep the hopping small, we need that the sepa- 
ration between vortices at the magnetic fields of interest 
is much larger than the coherence length. We therefore 
require that the field Hcz at which levels start to fill 
satisfies Hcz ^ Hc2- However, the temperature scale as- 
sociated with the field Hcz should not be too small since 
the physics of interest occurs only below that tempera- 
ture, suggesting materials with a relatively high Tc. A 
possible family of candidates that meet these criteria are 
the borocarbides . Another possibility is to work with 
a material that has H^z < Hd- Then even at magnetic 
fields just above iJci, when the vortices are very well iso- 
lated, a degenerate quasiparticle liquid will be formed in 
the vortex cores. Elemental Niobium, that has a rela- 
tively large Hd = 0.14T is possibly such a system. 

Throughout this work we have assumed that the su- 
perconductor is perfectly clean, but in any real system 
disorder is always present. The stability of the Luttinger 
liquid to weak disorder depends on the value of the inter- 
action strength g It is well-known that there exists 
a critical < Q such that for g > gc, any disorder kills 
the Luttinger liquid behaviour leading instead to a phase 
with localized quasiparticle excitations. Still, for suffi- 
ciently clean systems there is a range of temperatures for 
which the properties of the system are controlled by the 
LL fixed point, even though the ultimate zero tempera- 
ture state may be localised. In that case, the scanning 



tunneling conductance that we predict for the LL should 
be obtained in this crossover regime. 

On the other hand, if g < g^ the Luttinger liquid is, 
in fact, stable to weak disorder ||2^. Since g for the vor- 
tex core Luttinger liquid could be negative, the following 
interesting effect can occur. By varying the external mag- 
netic field, g can be varied and the system can be tuned 
through a one dimensional delocalization transition from 
the phase with localized quasiparticles (at small negative 
g) to the Luttinger liquid (at large negative g). 

The formation of the degenerate quasiparticle gas in 
the vortex core can be probed by measurements of the 
low temperature specific heat - a linear temperature de- 
pendence with a field dependent coefficient should obtain 
at the lowest temperatures. Once signatures of the for- 
mation of this quasiparticle gas are observed and some of 
its physical parameters (eg. the minigap and dispersion 
with kz) measured, it would be of great interest to look 
for possible Luttinger liquid behaviour in, for instance, 
STM tunneling into the vortex core. 

We also considered coupling the vortex quasiparticles 
to the collective modes of the vortex. We find that a 
vortex analog of the Peierls transition with a non-zero 
transition temperature could arise if the inter-vortex in- 
teractions are strong enough. 
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